Coulomb Drag for Strongly Localized Electrons: Pumping Mechanism 
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The mutual influence of two layers with strongly localized electrons is exercised though the 
random Coulomb shifts of site energies in one layer caused by electron hops in the other layer. We 
trace how these shifts give rise to a voltage drop in the passive layer, when a current is passed 
through the active layer. We find that the microscopic origin of drag lies in the time correlations of 
the occupation numbers of the sites involved in a hop. These correlations are neglected within the 
conventional Miller-Abrahams scheme for calculating the hopping resistance. 



Since the subject was first introduced in pioneering pa- 
pers , Coulomb drag between two parallel electronic 
layers has commanded a lot of attention. Stimulated by 
the early experiments the development of the the- 

ory proceeded in two directions: 

(i) progress on the formalism for calculating drag between 
conventional two-dimensional electron gases f^-jl^ . 
(a) accomodation of various realizations of interacting 
layers. These include electron- hole layers layers 
in the superconducting state electronic layers with 
tunneling links , strongly disordered layers at the on- 
set of Anderson localization , diffusive layers with cor- 
related disorder and double-layer systems in a per- 
pendicular magnetic field. In the latter case, the picture 
of drag depends on the magnetic-field regime, namely 
classically strong fields quantizing fields ||l9^^, the 
vicinity of the integer quantum Hall transition |21|, and 
the fractional quantum Hall regime p^-pst. 

For all these interacting two-dimensional systems, the 
theories of Coulomb drag shared the common scenario 
of quasiparticles (electrons, holes, or composite fermions 
p2|-p5t) in two contacting layers scattering off one an- 
other in the course of ballistic motion, diffusion or 
anomalous diffusion This scattering results in a 

non-zero average momentum transfer between the active 
(current-carrying) and passive (open-circuit) layers. A 
voltage drop is then induced across the passive layer to 
ensure the absense of a net momentum. 

Consider now the deeply insulating regime, where the 
localization radius of electronic states is smaller than the 
inter-electronic distance. Obviously, the momentum is 
not a good quantum number in this case, so that the 
conventional scenario of Coulomb drag does not apply. 
In addition, the picture of long-range time fluctuations 
of the electron density, and thus the language of spatially 
averaged vespoTise functions, is inadequate in the strongly 
localized regime. This is because the electron motion is 
due to hopping, which is characterized by an exponen- 
tially wide spread in the local transition frequencies. 

Instead of momentum exchange due to collisions, the 
coupling between the contacting strongly localized sys- 
tems is based on random shifts of energy levels in the 
passive layer caused by electron hops in the active layer, 



and vice versa. In this situation, it is not immediately 
obvious how these shifts "communicate" the overall di- 
rection of current from the active to the passive layer. 
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FIG. 1. Schematic drawing of (a) active layer with applied 
voltage V and two localized states 1 and 2. U indicates the 
Coulomb energy by which level 1 (2) is elevated if level 2 (1) 
is occupied, (b) The passive layer also involves two localized 
states 1' and 2' , but has no voltage applied. 

This question is addressed in the present paper. We 
find that Coulomb drag in the strongly localized regime 
is due to the correlated character of hopping transport. 
To clarify this point, consider two neighboring sites, 
1 and 2, which belong to the current-carrying cluster 
p6| . At nonzero temperature the time evolution, n\(€) 
and n2{t), of the occupation numbers of sites has the 
form of telegraph noise. In equilibrium, the average 
Ki2{t) — {ni{t)n2{t -l- r)) is an even function of r. Sup- 
pose now that current flows in the direction 1 — > 2. Then 
we have Ki2{t) > Ki2{—t) for r > 0. This reflects the 
fact that, as the hops within a pair of sites occur pref- 
erentially from 1 to 2, the occupation numbers change 
in a certain sequence. This asymmetry is the analog of 
the current-induced asymmetry between wave vectors q 
and —q of the thermal density fluctuations in the metal- 
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lie regime, and thus, it is responsible for Coulomb drag 
between layers with strongly localized electrons. 

Note in passing, that conventional theories of hopping 
transport neglect the asymmetry in Ki2{t). In both 
the non-interacting and the interacting ||2^ cases, 
the resistance of an elementary hop is computed under 
the assumption of uncorrelated occupation numbers, i.e. 
{ni{t)n2{t')) = (m) (na). 
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FIG. 2. Schematic time evolution of occupations ni{t) and 
722 (t) of levels 1 and 2, respectively, in the active layer under 
conditions of a finite current. 

To show how the preferential sequence in the change 
of the occupation numbers transforms into drag we con- 
sider the simplest possible model as illustrated in Fig. I. 
Within this model, the active and the passive layer are 
each represented by a pair of sites (1 and 2) coupled to 
metallic contacts {I and r). Both "active" and "passive" 
pairs of sites can be either empty or singly occupied. The 
corresponding conditions are £1^2 + U > Ep + eV for the 
active layer, where U is the Coulomb interaction between 
two electrons on sites 1 and 2 and V the applied voltage, 
and ei',2'+U > Ep for the passive layer. Due to this con- 
dition an electron occupying, say, site 2 elevates the en- 
ergy level of site 1 above the Fermi level, thus forbidding 
the tunneling process I 1 (see Fig. 1). We will treat 
this model first in the strongly-nonlinear regime, in which 
the temperature T is much lower than V . Subsequently, 
we will consider the Ohmic regime, where 1^ <C T. 

In the strongly-nonlinear regime, we can neglect all ac- 
tivation processes and thus, only a single sequence of hops 
is possible in the active layer. Within this sequence, the 
occupation numbers ni and n2 of sites 1 and 2 undergo 
the transformations (1,0) (0,1) (0,0) (1,0). 
With each repeated cycle of this sequence, an electron is 
transferred from the left to the right contact. The time 
evolution of ni and 712 during one cycle is illustrated in 
Fig. 2. The strong asymmetry in Ki2{t) is evident. The 
average current through the active layer is equal to [ p8|j29| 
{la) = s./ (ji + T2 + T^), where ti, T2 and T3 are the av- 
erage waiting times for the transitions Z — * 1, 1 — > 2 and 
2 — > r, respectively. In the strongly nonlinear regime, the 
current {Ig) is independent of the voltage drop, V. 
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FIG. 3. Variants of evolution of the level occupations in 
the passive layer during a single cycle in the active layer. 

Consider now the response of the passive layer to a 
single cycle of (ni(t), n2{t)). In the configuration ni ~ I, 
n2 — 0, the occupied site 1 elevates the level £i' by some 
Wi , and the level £2' by some W2 = Wi - A < Wi. Con- 
versely, when ni = 0, 71-2 = 1, the level ey is elevated by 
W2, while the level £2' is elevated by Wi. (We define the 
level positions £i', £2' in the passive layer with respect to 
an "empty" active layer, ni = n2 — 0.) As we will see, 
the pair 1', 2' is most "sensitive" when the conditions 



Ep > £1/ +Wi> £2' -t- W2 
£2' +Wi> Ep 



(1) 
(2) 



are met. Indeed, by virtue of Eq. (^, the passive layer 
can only make the transition 1' — > 2' during the interval 
ti < t < t2 when site 1 in the active layer is occupied. 
This is illustrated in Fig. 3a. If the hop 1' 2' does 
take place, then two transitions within the passive layer 
become possible during the next time interval t2 < t < t^: 

1) the electron on site 2' goes back to 1', cf. Fig. 3b. This 
transition is energetically favorable due to Eq. (^. 

2) the electron tunnels from 2' into the right contact, r', 
emptying the passive layer. Since ei> + W2 < Ep, this 
opens the possibility for the hop /' 1' (see Fig. 3c). 

In the first case, there is no transfer of charge from the 
left to the right contact in the passive layer during the 
cycle in the active layer. In the second case, however, 
provided that both transitions 2' r' and V — > 1' took 
place during the interval t2 < ^ < ^3, the cycle in the 
active layer results in the transfer of an electron from V 
to r' . This transfer is nothing but the drag current. 
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The way the drag current is induced here closely re- 
sembles the operation of a classical electron pump \ 30 - 3^] 
where two phase-shifted rf signals are applied either to 
the barriers or to the gate electrodes ||3l|,^ of a sin- 
gle or multiple quantum dot structure. In contrast to 
adiabatic quantum pumping |3^,Q , classical pumping is 
due to the Coulomb blockade, which forces an electron 
from one contact to enter the dot, and then to leave the 
dot into the other contact during each rf cycle. In the 
mechanism of hopping drag considered above, the pulses 
ni{t) and 77-2 (<) play the role of the rf signals, their mu- 
tual phase shift being governed by the direction of the 
current in the active layer. 

We now turn to the calculation of the average drag cur- 
rent, (Id)- Denote with Ti,72, and the waiting times 
corresponding to the transitions 1', 1' — > 2', and 

2' —^ r' in the passive layer. The simple system consid- 
ered here mimics drag between two-dimensional hopping 
layers when the transitions 1^2 and 1' 2' consti- 
tute "bottlenecks" for the transport in active and passive 
layer, respectively, i.e. T2 and T2 are the "long" waiting 
times. In this limit we have (/a) « e/T2. The calcula- 
tion of (Id) is greatly simplified if the "short" waiting 
times are related as follows 7i <C ^ T3. The latter 
conditions ensure that the passive layer returns to the 
ground state after each cycle in the active layer. Indeed, 
if the transition 1' 2' took place during the time in- 
terval ti < t < t2 (the corresponding probability is equal 
to pi'2> = 1 — exp[— (^2 — ii)/^], then the probability 
that, during the subsequent interval ^2 < ^ < ^3, both 



transitions 2' 



and I' — > 1' take place is close to 1. 



This is because the characteristic duration of this inter- 
val (^3 — ^2) ^ T3 is much longer than and 7i. In 
principle, the transition 2' r' opens the possibility for 
two follow-up hops, namely, I' — > 1' and r' — > 2'. How- 
ever, since 7i ^ 73, tunneling of an electron from the left 
contact onto site 1', thereby returning the passive layer 
into its ground state, occurs before the back hop r' ^ 2'. 
Once the transition I' V occurred, the passive layer 
remains in the ground state until the end of the cycle in 
the active layer t = t4 (see Fig. 2). Thus, calculating the 
drag current reduces to averaging pi>2' which yields 



T2 



T2 + T2 



{Ia) = 



T2 + T2 



(3) 



Remarkably, the drag current is not exponentially small 
compared to the current in the active layer. Instead, 
both currents are of comparable magnitude. This result 
is clearly very different from naive predictions based on 
spatially-averaged response functions such as a widely 
used Fermi-golden-rule-type expression for the drag re- 
sisitivity 0. 

In principle, hops in the passive layer might, in turn, 
affect transport in the active layer. In particular, occupa- 
tion of site 2' shifts level 2 upward. As a result, the hop 
2 — » 1, in the direction opposite to the net current, might 



become energetically favorable. The condition that such 
"feedback" does not occur is ei > £2 + A. As a conclud- 
ing remark on our simplified model, we note that for a 
given polarity of voltage across the active layer (see Fig. 
1) the activationless drag exists only if e'^ < e'2. 

We now turn to the linear regime V where trans- 
port in the active layer becomes more complicated in two 
respects. Firstly, charge-transfer processes become possi- 
ble in both directions, I — *■ r and r I, with only a small 
difference in their frequencies due to the applied voltage. 
Secondly, the dynamics involves "round-trip" processes 
such as^— >l^ZorZ^1^2^1^Z that do not 
result in charge transfer between the contacts. 

It is important to note that the traditional description 
of hopping transport ||2^] , based on the Miller- Abrahams 
network, does not capture the realistic occupation dy- 
namics of the sites 1 and 2. Indeed, in this description, 
a pair of sites is considered as an isolated system so that 
the characteristic times for the changes in the occupa- 
tions ni and 712 are long (~ T2 in our notations). This 
is indeed the case for the non-ohmic regime considered 
above where an electron has to spend a time T2 at site 
1 before the transition 1 — > 2 takes place, since it cannot 
return to the contact I without activation. In the ohmic 
regime, however, while the waiting time, T2, is long, the 
occupation of sites 1 and 2, constituting a "bottleneck", 
changes many times between successive hops 1 — > 2 (or 
2^1). In other words, a typical charge-transfer process 
is preceded by many round-trip processes. 

As in the non-ohmic regime, we restrict our considera- 
tions to a particular domain of parameters for which the 
analysis of drag is greatly simplified. We assume that (1) 
the energies of both sites in the active layer are above the 
Fermi level (^1,2 > Ep). Thus their average occupations 
are small, {ni{t)) , {n2{t)) <C 1, i.e. the states 1 and 2 are 
empty for long periods before an electron enters for a time 
of the order of ri or T2. In addition, we assume for the 
passive layer that (2) the sites 1' and 2' are symmetric in 
space, so that 7i « 73, and close in energy (within tem- 
perature). As a result, unlike the situation in Fig. 3a, the 
site 1' is also elevated above the Fermi level, when the site 

1 in the active layer is occupied. Finally, we demand that 
(3) the "long" time T2 in the passive layer is shorter than 
the typical time interval during which the active layer is 
empty, but longer than the time when one of sites 1 or 

2 are occupied. The latter assumption implies that an 
electron in the passive layer hops many times from 1' to 
2' and back between successive round-trip processes in 
the active layer. We also demand 71,72 ^ ti,T2. 

Under these conditions, we find a remarkable result, 
namely that the drag current is simply half the current 
in the active layer. To understand how this comes about, 
consider a single charge-transfer process in the active 
layer, Z — s- 1 — s- 2 — > r. Imagine first that the passive 
layer is in the state ni' = 0, n2' = 1 when the transi- 
tion I — > 1 takes place. In this case, nothing happens in 
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the passive layer until the transition 1 — > 2 takes place, 
after which the passive layer will undergo the transition 
2' r' in essentially all cases. This will be immediately 
followed by V ^ V . With the transition 2 — > r the system 
effectively returns to its initial state and one electron has 
been transferred in both layers. Since the passive layer 
has a 50% chance of being in the state ny = 0, n2' = 1 
at the outset, this implies 

{I<l)-l{Ia). (4) 

Eq. (Q) is based on the fact that there is no comparable 
drag current when the passive layer initially has an elec- 
tron in 1'. In this case, I — > 1 causes the passive-layer 
transitions 1' — > /' and then r' —>■ 2' , which are reversed 
after 1^2 occurs in the active layer. Using similar rea- 
soning, it can be verified that round-trip processes in the 
active layer do not lead to a net drag current. 

An actual sample in the deeply localized regime will 
consist of a network of transresistors of the type described 
above. Hence, we should compute the transresistance of 
this network. This is by no means a trivial task, even 
if the transrcsistances between elements of the conduct- 
ing networks in active and passive layer are known. This 
can be illustrated by considering the way in which tran- 
srcsistances combine when connected in sequence or in 
parallel. These situations are readily analyzed in terms 
of Kirchhoff's laws. Two transresistors in sequence have 
a transresistance equal to the sum of the individual trans- 
resistances, 

R^'^^R^'+R?' (5) 

similar to ordinary resistors. Two transresistances in par- 
allel are less simple. Here, one can show that 

(1||2) ^ RiR'iRI^^ + R2R'2R^t^^ 

* {Ri+R2){R'^+R'2) ' ^ ^ 

where Ri {R[) denote the resistances of the active (pas- 
sive) layer of the i-th transresistor. These expressions al- 
ready allow one to draw an important conclusion about 
the net transresistance of the entire network. Even 
though the resistances Ri and R[ are exponentially large, 
adding transresistances in sequence or in parallel does not 
lead to exponential changes in the transresistance. In this 
sense, the drag current is of comparable magnitude as the 
current in the active layer, even for the entire network. 

In conclusion, the analysis of particular realizations of 
the active and passive layers, carried out in the present 
paper, illuminates the physics underlying the strong drag 
in the localized regime. The discreteness of the hopping 
electrons gives rise to Coulomb blockade, which, in turn, 
opens the possibility of classical pumping. 
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